Nonequilibrium Dynamics Across an Impurity Quantum Critical Point 
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Whether a small quantum mechanical system is able to equilibrate with its environment once an 
external local perturbation drives it out of thermal equilibrium is a central question which cuts across 
many different fields of science. Here we consider such a problem for a correlated quantum impurity 
coupled to a fermionic reservoir and driven out of equilibrium by local quantum quenches such as 
those recently realized in optical absorption experiments on single quantum dots. We argue that 
equilibration in this problem is deeply connected to the occurrence of Kondo Effect at low energy 
and that a highly non trivial dynamical behaviour may emerge whenever a local quantum critical 
point intrudes between a conventional Kondo screened phase and a Kondo unscreened one. We 
discuss this issue in the context of the Anderson Impurity model coupled to a pseudo-gap reservoir 
by using a correlated time dependent variational wave function that is able to qualitatively describe 
this physics 
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Introduction - Small interacting quantum mechanical 
systems coupled to an external environment represent 
basic paradigms of transport, dissipation and non equi- 
librium phenomena. Understanding their dynamical be- 
haviour and to what extent those system are able to equi- 
librate with their environment when a local perturbation 
drives them far from equilibrium is therefore crucial in 
many different physical contexts. While for systems cou- 
pled to a bosonic environment this question has been 
around since the early days of Caldcira-Leggctt models 1,2 
and it is keep attracting a lot of interest*^—, much less is 
generally known about the out of equilibrium relaxation 
dynamics of small quantum systems coupled to fermionic 
reservoirs. Experimental realizations of these fermionic 
quantum impurity models traditionally involve diluted 
magnetic impurities in metals^, quantum dots and single 
molecules attached to leads^— or magnetic adatoms on 
metallic surfaces 12 . While well suited to study transport, 
these setup do not probe directly the nonequilibrium 
relaxation dynamics. Recent achievements in control- 
ling optical properties of single semiconducting quantum 
dots brought a new twist in quantum impurity physics. 
By shining coherent light on those systems, experiments 
are able to probe the out of equilibrium dynamics in- 
duced by a sudden switching of a local perturbation, a 
so called local quantum quench^— . The experimental 
evidence of the Kondo effect, a smoking gun of strong 
local correlations, in the absorption spectrum of a sin- 
gle self assembled quantum do1ji£ opens the exciting pos- 
sibility to explore in a highly controlled setup the non 
equilibrium dynamics of quantum impurity models cou- 
pled to fermionic environments. A fundamental question 
this experimental breakthrough brings over is whether 
and under which condition a quantum impurity is able 
to equilibrate with its fermionic environment after a local 
perturbation, such as a sudden switch-on of the coupling 
with the bath, drives it far from equilibrium. A very well 
studied example along this line is the Anderson Impu- 



rity Model (AIM)i^ a paradigmatic model describing a 
single interacting electronic level coupled to a fermionic 
reservoir. A number of theoretical investigations on its 
out of equilibrium dynamics found that quite generically 
thermalization after a local quantum quench occurs22r— 
. Yet, when the system is deep in the Kondo regime an 
exponentially long time scale controls this relaxation, as 
the system flows toward the strong coupling (SC) fixed 
point. The Kondo effect is, however, not the only possi- 
ble fate for a quantum impurity at low energy. A different 
and perhaps more intriguing possibility, which motivates 
the present paper, may arise whenever a mechanism com- 
peting with the Kondo effect is present, which is able to 
drive the system away from the SC fixed point and even- 
tually to renormalize to zero the effective coupling with 
the bath. This is what happens when an impurity quan- 
tum critical point intrudes between a SC Kondo screened 
phase and a Local Moment (LM) phase. Here we study 
the out of equilibrium real time dynamics across such an 
impurity quantum criticality. We consider the simplest 
example showing this competition, namely the Anderson 
impurity model coupled to a pseudogap reservoir— ~—. 
By using a time dependent correlated variational wave 
function we show that an highly non trivial dynamics 
emerges as a result of the complex low energy flow of the 
model. 

The Model, the Local Quantum Critical Point and its 
variational description - The model we are going to con- 
sider is the so called pseudogap Anderson Impurity Model 
(pg-AIM) describing a single spin-full fermionic level cou- 
pled to a reservoir of non interacting electrons 

HAIM — Hbath + Himp + Hhyb = ^ £k /k CT /ktr + 

ko- 

+ Y("- 1 ) 2 + E y k (4/k<r + fc-c) (1) 
k<r 

with n = c\c a . The properties of the bath 
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are encoded in the hybridization function T(e) — 
n ^2 k 5 (e — Ek) , that we assume to vanish as a 
power law at the Fermi energy e = 0, namely T (e) — 
(1 + r) r \e/W\ r 9(W - \e\). The equilibrium phase di- 
agram of the pg-AIM is extremely rich (see Ref i i 28 
for reviews), featuring at particle-hole symmetry and for 
0<r<l/2,a quantum phase transition at a criti- 
cal value of interaction U c between a Kondo screened 
phase (for U < U c ) and a Local Moment (LM) regime 
(for U > U c ) where the impurity becomes asymptoti- 
cally free at low energies. As opposite, for r > 1/2 the 
only stable fixed point is the local moment one and no 
Kondo effect can be stabilized for an Anderson Impurity 
in such a gapless reservoi r 25 i 26 i 29 . The equilibrium prop- 
erties of the Kondo to Local Moment quantum criticality 
have been the subject of a large literature (see Ref. 30 
for a recent work). Quite remarkably though, its out of 
equilibrium dynamics is an almost unexplored field 3 ^. 

As it will be relevant for our study of the dynamics we 
start our discussion by presenting a simple and intuitive 
equilibrium variational description of the half filled pg- 
AIM, in the spirit of the Gutzwiller ansatz, that is able 
to qualitatively describe salient features of this quan- 
tum phase transition. We write this correlated ansatz 
as |\E r ) = V\$z) where \§>z) is the groundstate of an 
half-filled non interacting (i.e. U = 0) resonant level 
model (RLM) with a variational hybridization y/ZVk 
while V acts on the impurity Hilbert space and reads 
V = J2n=o 12 | n}(n\. Here | n) { n| projects onto con- 
figuration at fixed impurity charge n, while A„ are varia- 
tional parameters. The key idea behind the above ansatz 
is that (i) the coupling between the impurity and the 
reservoir is the crucial quantity to distinguish between 
the two sides of the QCP and (ii) strong local correlations 
can substantially renormalize this effective coupling. In 
order to solve the variational problem we should minimize 
the energy functional E = (\P | H \^>) with respect to all 
the variational parameters . At half filling one can show-2 2 - 
that only one independent variational parameter enters, 
that we choose to be \[Z. This is the so called quasi- 
particle renormalization factor, which in equilibrium also 
controls the weight of the low-energy peak in the spectral 
function. The saddle point equation reads 

^= + ^J d smeA z{ s) = 0, (2) 

where Az (e) is the spectral function of a RLM with a 
renormalized hybridization \fZVk ■ Its solution as a func- 
tion of U /r and for different bath exponents r is shown in 
figure [TJ When r = (metallic bath) Z smoothly crosses 
over to an exponentially small value in the strong cou- 
pling Kondo regime U T, Z ~ exp(— it U/16T). How- 
ever as soon as r ^ it is possible to have a vanishing 
quasiparticle residue at finite U. The critical interaction 
strength U c for this to happen can be obtained in closed 
form U c = MElAir and correctly disappears as r — > 0. 
For U < U c the variational ground state has a finite Z, 
meaning the system is a local Fermi Liquid as one expect 
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FIG. 1: Gutzwiller results for the pg-AIM at half filling. We 
plot the renormalization factor Z as a function of U /Y for 
different values of r, the power-law exponent of the bath. We 
see that upon increasing r from the metallic case a critical 
point emerges which separate a Kondo SC phase (Z 7^ 0) 
from a LM (Z = 0) fixed point. Inset: Critical interaction U c 
as a function of r. 



at the Kondo strong-coupling fixed point. As opposite 
for U > U c the impurity is effectively decoupled from 
the bath, since as Z — !• also the effective hybridization 
vanishes. This is a sign the system flows to the LM fixed 
point. Close to U c the quasiparticle weight vanishes lin- 
early, Z ~ U c — U. While the above variational analysis 
provides a simple picture for the quantum phase transi- 
tion between a Kondo (Z ^ 0) and a LM (Z = 0) regime, 
a more careful treatment of quantum fluctuations would 
have revealed an even more rich phase diagram. Indeed 
it has been shown by means of NRG2£ and field theo- 
retical analysis^ that such a QCP is present only for 
r £ (0, 1/2), while for r > 1/2 the system always flow to 
the local moment fixed point. This is not captured by 
our simple variational wave function which always shows 
a finite U c for all r > 0. As a result we will confine our 
analysis to the regime < r < 1/2. 

Dynamics after a local quantum quench - Let us now 
consider the dynamics induced in the pgAIM by a lo- 
cal quantum quench. In particular we will focus on two 
classes, namely a sudden switching of (i) the hybridiza- 
tion 14 or (ii) the interaction U. Since in a local quantum 
quench the energy injected into the system does not scale 
with the size we expect that in the thermodynamic limit, 
if the system thermalizes, the dynamics will be relaxing 
toward the ground state of Haim^- In order to study 
the role of QCP on the dynamics of pgAIM we resort here 
to a time dependent strongly correlated Gutzwiller vari- 
ational approach, following recent developments ^ 35 ' 37 . 
To this extent we introduce an ansatz for time-dependent 
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FIG. 2: Variational dynamics for the pgAIM with U/T = 9, 
r = 0.4 after an hybridization (full lines) or interaction 
(dashed lines) quench. We plot the impurity double occu- 
pancy D(t) and its conjugate variable tp(t) (see inset), the 
renormalized quasiparticle weight Z(t) and the hybridization 
energy Eh(t). We see that independently from the initial con- 
dition the dynamics relaxes to the variational groundstate of 
the pgAIM model (as shown by the arrows), namely the sys- 
tem thermalizes after a local quench. In the bottom right 
panel we plot D(t) for different values of U/V = 5,7,9,11 
from top to bottom. 



FIG. 3: Variational dynamics for the pgAIM (r = 0.4) after a 
local quench across the SC-LM quantum critical point. Upper 
panel: dynamics of double occupancy D(t) and quasiparticle 
weight Z(t) for an interaction (red) or hybridization (black) 
quench, with U/V = 25. In both cases local quantities do not 
relax. Bottom panels: the average hybridization between the 
impurity and the bath shows coherent undamped oscillations 
(left) for both type of quenches. This can be traced back 
to the vanishing at U c of the energy scale rj controlling the 
coupling between the impurity sector and the bath modes (see 
right and main text). 



many body wave function 

Mt)) = v(t)Mt)), (3) 

where now V{t) = J2 a X n (t) e llfn ^ \ n)(n\ is a time de- 
pendent operator acting on the impurity Hilbert space 
with X n ,ip n being variational parameters, while |$(£)) is 
the solution of a time dependent Schroedinger equation 
with renormalized non interacting hamiltonian whose 
form will be specified in the following. The variational 
dynamics is then obtained by making the following real 
time action 

S = J dt(*(t)\id t -H\*(t)) (4) 

stationary with respect to all the variational parameters 
X n ,ip n , | $(£)). We quote the final result of this calcula- 
tion, which goes along the same line of previous case^. 
Assuming to start from a particle-hole initial prepara- 
tion for the dot and the reservoir, then the dynamics for 
double occupancy D(t) and its conjugate variable (f(t) 
reads 

<p = 4E hyb (t) vAD (1/2-D) sirup (5) 
D = U + E hA t) ^ 2 D _ D) cos V (6) 



with E hyb {t) = {<b(t)\H hyb \$(t)), while the state |$(t)) 
satisfies the Schroedinger equation of a RLM 

id t |$(t)) = (H bath + ^Z(J)H h y b ) I $(*)> (7) 

with a time-dependent hybridization W Z(t) — 
4 y/D{t) (1/2 — D(t)) costp(t) . It is important to stress 
at this point that the above dynamics fully accounts 
for the coupling between the variational parameters 
D, ip entering the projector and the fermionic degrees 
of freedom in the wave function Let us now 

discuss the solution of the variational dynamics (JSJ?]) by 
firstly considering local quenches that end in the Kondo 
screened phase, namely for U < U c . In figure [5] we plot, 
for the two kind of local quenches we are interested 
and for a bath exponent r = 0.4, the dynamics of three 
relevant quantities namely the double occupancy D(t), 
its conjugate phase ip(t) and the hybridization energy 
Eh yb (t). We see that after an initial transient, for time 
scales larger than some relaxation time r re i the dynamics 
reaches a steady state which is independent from the 
initial condition. A direct comparison shows that such a 
steady state coincide with the (variational) ground state 
of the pseudo-gap AIM, that we plot as a dashed line in 
the figure. In other words the system thermalizes after 
a local quench. However, a very different dynamical 
behaviour emerges when the system is quenched deep in 
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the LM regime, as we show in figure [3j In particular the 
double occupancy shows very fast oscillations around 
an average value which is different from the expected 
ground state value. More interestingly, the hybridization 
strength, which measure the effective coupling between 
the impurity and the bath, starts oscillating around 
zero. This enlights the fact that above the critical point, 
as the impurity effectively decouples from the bath, the 
dynamics cannot equilibrate. 

Discussion - The thermalization of the pgAIM for 
U < U c is consistent with the behaviour of the conven- 
tional AIM in a metallic bath 2 ^— . It reflects the fact 
that, as the low-energy /long-time limit is approached, 
the model flows toward the SC fixed point where the 
impurity stays always coupled to the reservoir and the 
system behaves as a local Fermi Liquid- . This phase 
is adiabatically connected to the weak coupling RLM 
whose dynamics after a local quench leads to thermal- 
ization. A closer look to the variational dynamics gives 
further insights into this thermalization process. Indeed 
one can explicitly show 3 -^ that the equilibrium ground- 
state variational solution is a steady state (SS) of the 
dynamics (|5I6I7I) . This is immediately seen for the classi- 
cal degrees of freedom by noticing that (i) (p = implies 
sinip* = in the SS and (ii) after simple algebra the 
equation © becomes the saddle point condition for 
Z+, provided the hybridization energy of the RLM ther- 
malizes at long time, Eh y b{t ^> r re i) — E* lyb , as one can 
show to be the case 3 -^ . This result suggests to linearize 
the dynamics around the SS and to extract information 
on the relaxation time upon approaching the QCP. A 
simple but tedious calculation gives the following reduced 
dynamics^ for the double occupation 5D(t) = D(t) — D+ 

SD = a5D + r] SE hyb (8) 

in terms of 5Eh y b = Ei iy b(t) — E* lyb , the deviation from 
steady state of the average hybridization. The coeffi- 
cients a and rj give information on characteristic relax- 
ation timescales and are purely equilibrium steady state 
properties. They read respectively a ~ E i ^ jh /Z i , and 
7] ~ (1 — 4D*) E* lyb . Their behaviour with U, as plotted 
in figure [3J reveal that while a is finite (and negative) 
approaching the QCP as E* lyb ~ \fZl, the coupling r] 
between the impurity sector and the bath modes van- 
ishes at U c as 77 ~ ^fZl ~ VU C — U. This vanishing 
energy scales suggests a diverging time-scale for thermal- 
ization (and relaxation) of local observables. As a result 
for quenches above the QCP no steady state can be ob- 
tained within the above variational dynamics, as already 
revealed by the numerical solution plotted in figure [3J 

It is interesting to discuss whether the above picture 
is robust against quantum fluctuations that are not in- 
cluded in the Gutzwiller variational scheme. We expect 
these to play a special role on the LM side of the QCP, 
to damp out the oscillations and eventually lead to a 
(non-thermal) steady state. In addition we expect that 
around U c , and particularly for quenches ending at the 



QCP, important effects beyond mean field will be crucial 
to capture correctly the quantum critical relaxation dy- 
namics. A possible framework to treat non-equilibrium 
effects beyond this time dependent variational scheme is 
offered by the slave spin formulation of the AI M 35 ' 39 , 
where the impurity electron is represented as an auxil- 
iary fermion coupled to a two-level system (TLS). Quite 
interestingly , the Kondo effect emerges in this picture 
as a localized phase of an effective sub-ohmic spin bo- 
son model^ 9 -. As far as the non equilibrium dynamics is 
concerned, while a simple mean field decoupling of the 
TLS and fermionic sectors leads to an evolution which 
is fully equivalent to the time-dependent Gutzwiller dis- 
cussed in this paper, quantum fluctuations beyond the 
saddle point can be included systematically. These re- 
sult into coupled dynamical equations for the auxiliary 
fermionic and TLS self-energies, similar to large N ap- 
proaches which are known to capture some aspects of 
the quantum critical relaxation dynamics 27,40 . These is- 
sues represent interesting directions for further investiga- 
tions. In spite of that, we could say that the qualitative 
picture of the non equilibrium dynamics, including a di- 
verging time scale and the absence of thermalization for 
quenches in the LM regime appear as a robust features 
of the pgAIM model that are well captured by our simple 
variational wave function. In this respect it is interesting 
to notice that a recent investigation of the nonequilibrium 
dynamics in the Ferromagnetic (FM) Kondo mode l 41 ' 42 , 
that always flows to the LM fixed point, finds a similar 
lack of ergodicity in the long time dynamics. 

Finally, we conclude by discussing possible connections 
with experiments. As we mentioned in the introduc- 
tion, recently the non-equilibrium dynamics of an AIM 
after a local quench has been addressed, indirectly, in 
experiments on single self- assembled quantum dots 1 ^ by 
measuring the absorption/emission spectrum after opti- 
cal excitations. This spectrum turns to be related to 
the statistics of the work done during the quenc h 14 ' 15 , 
which features a clear signature of the Kondo effect in 
the form of an edge singularity at low values of the work. 
A natural question would be therefore whether one can 
arrange a system with a similar Kondo to Local Moment 
QCP in these systems and how to characterize this quan- 
tum phase transition from the point of view of the non- 
equilibrium dynamics after a local quench. As far as the 
first point is concerned, one obvious possibility would be 
to play with the conduction electrons density of states, 
so to arrange for a semi-metallic behaviour. Since the 
conduction electrons are provided by a doped semicon- 
ductor (typically n-doped GaAs), this may be hard to 
realize in practice. A more interesting scenario would be 
to probe the dynamics of a system of two Kondo spins 
anti-ferromagnetically coupled, which is known to posses 
a similar quantum phase transitions 3 -. In that case, a 
clear dynamical signature of the QCP would emerge by 
looking at the asymptotic behaviour of the optical spec- 
trum, where the X-ray edge singularity is expected to be 
cut-off for a much weaker (logarithmic) singularity. 
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Conclusions - In this paper we have discussed the 
nonequilibrium dynamics after a local quench in the 
pgAIM, as a prototype of fermionic impurity quantum 
critical point. We have argued that equilibration in this 
problem is deeply connected to the occurrence of Kondo 
Effect at low energies. While for local quenches in the 
Kondo phase thermalization occurs, this is not the case 
when quenching across the QCP deep in the LM regime. 



This is correctly captured by a simple and intuitive time 
dependent variational wave function which in addition is 
able to describe the divergence of the relaxation time as 
the QCP is approached. 
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A. Gutzwiller wave function for the AIM 

Here we briefly sketch the main results of the 
Gutzwiller variational approach to the AIMi~— . We write 
the hamiltonian as 

HaIM = Hbath + Him P + Hhyb = ^ f^ a /k<x + 

+f (n - if + J2 V * (4 A. + h.c.) (9) 

kcr 

with n = <V a drx and consider a wave function of the 
form 



m = v\$> 



(10) 



where is the groundstate of a non interacting (i.e. 
U = 0) half filled resonant level model (RLM) with a 
variational hybridization A Vk while V acts on the im- 
purity Hilbert space and reads V = J2n=o l 2 ^« I n )( n \- 
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Here | n){ n\ projects onto configuration at fixed impurity 
charge n, while A„ are variational parameters, related to 
the correlated probability of having a given charge n on 
the impurity. This latter quantity is defined as 



Pn 



(tf|ra)H*) = \ 2 n ($ X \n)(n\$ X ) 



(11) 



We can impose without loss of generality the normaliza- 
tion condition on the wave function 



(12) 



that also implies ^ n P n = 1. In addition one can easily 
convince that for the above wave function the following 
relation holds at particle-hole (PH) symmetry 



(® x \v 2 m\) = ($aH*a) = i, 



(13) 



or equivalently n Pn = 1- Since at PH symmetry 
we have Pq = Pi = D from these two latter conditions 
we conlcude also that Pi = 1 — 2D. In order to solve 
the variational problem we should minimize the energy 
functional E var = | H |\&) with respect to all the vari- 
ational parameters. The identity (|13[) greatly simplifies 
the evaluation of the variational energy. Indeed using 
this result one can show that 4 



bath\®\ 



as well as 



(*\H hvb \*)=VZ($ x \H hyb \$ 



(14) 



(15) 



where \J Z is a function of P n that at half filling reads 
^/Z = (1/2 — D). Using these results we can write 

the variational energy as 

E var = <tf = ($a | H bath \$ x ) + 

+VZ{$x\H hyb \$ x ) + UD (16) 

Variation with respect to |$a) at fixed D gives a renor- 
malized RLM hamiltonian with A = \[Z, The sad- 
dle point equation for Z is then obtained by taking 
the derivative of the variational energy. Using the 
fact that for a RLM with hybridization \~ZVk we have 
d Egg/dVz = (Hhy b ) gs and by writing this latter aver- 
age as 



ir^/Z 



def(s)eA z (e) (17) 



we obtain the equation in the main text, where we have 
defined Az (e) as the spectral function of a RLM with a 
renormalized hybridization \[ZV k 



Az (e) 



ZT(e) 



(e-ZA R (e)) 2 + (ZT(e)f 



(18) 



with 



„ f duj T(u>) v — • 
A fl (e) = P / — -±-L T(e) = TrJ2v k 2 5(e~e k ). 

J k 

(19) 



I. TIME DEPENDENT GUTZWILLER FOR 
THE AIM 

Extensions of the Gutzwiller wave function to the time 
dependent case have been recently introduced, both for 
lattice models^ and for quantum impurities 6 . The evalu- 
ation of the real-time action on the time dependent wave 
function 



|*(t)> = v{t) |$(t)), 



(20) 



where now V(t) = J2 a A n (t)e^»W \ n){n\ while 
is the time dependent wave function of a proper non in- 
teracting RLM, can be done along similar lines as in the 
equilibrium case. After imposing the normalization con- 
dition 



m)\r 2 (t)mt)) = i, 

and noticing that at PH simmetry we have 

(<£(t)\-P 2 (t)nMt)) = m)\nMt))=l 
we obtain the following real-time action 

S = J dt(V(i)\id t -H\y(t)) = 
= Srlm+ I dt (iipD-UD) 



(21) 



(22) 



(23) 



where Srlm is the action of an effective time dependent 
RLM that is 



Srlm 



dt($(t)\id t -H RLM (t)mt)). (24) 



with hamiltonian H R LM(t) — H b ath + \J Z(i) Hh yb where 
the time dependent renormalization factor \fZ~ reads 



Vz = 4y/D(l/2~D) 



costp 



(25) 



From the action (|23|) we get the variational dynamics 
for D{t) and tp(t) 



D = 4E hyb (t) ^/D (1/2 - D) sirup 



(26) 



1 - AD 

ip = U + E hyb (t) - == costp (27) 

v /D(l/2-£>) 

and the Schoedinger equation for 

idt |$(t)) = H RLM (t)\ $(t)) . (28) 

In order to solve the dynamics of the effective time de- 
pendent RLM we can proceed with an equation of motion 
appoach as in — . Let us define the following averages 



N, 



kp 



(29) 



whose dynamics is obtained by evaluating the commuta- 
tors with H RLM {t) 



: ,W k = -e k W k +J2v p N kp -V k 



k n 



(30) 
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* Nk P = (e P - £k) N kp + V p W k - V k W } 



in = 



(31) 
(32) 



where we have defined V k (t) = ^/Z{t)V k . We note that, 
since PH symmetry is preserved at any time even for a 
time dependent hybridization, the occupation of the dot 
is fixed and the last equation is irrelevant, asn = 1. The 
hybridization strength entering the classical equations of 
motion read 

Ehyb(t) =J2v k {W k + h.c.) (33) 



A. Long time dynamics and Thermalization 

We now want to show that the equilibrium groundstate 
variational solution is a fixed point of the above dynam- 
ics. In order to see this let us start from the classical 
dynamics and search for a stationary solution, D 
By putting derivatives to zero we get 



sirup* = 



and 



= U + E, 



1-4D* 



hyb 



vAD* (1/2 - AO 



where we have defined 



Ehyb 



lim. 



Ehybit) 



(34) 



(35) 



(36) 



as the steady state limit of the hybridization energy. 
Equation (|35p is nothing but the equilibrium saddle point 
energy condition for D 



dE9 s 
dD 



U 



BE 9S 

Urj hyb 

dD 
-- U 



u 



^hyb 



dE^bdVZ 
dVz dD 
I- AD 

y/D(l/2-D) " 



(37) 
(38) 



provided that the steady state value of the hybridization 
energy, E£ b , coincides with the corresponding ground 

state value Ef^ b as a function of \fZl, namely 



(39) 



Therefore in order to proceed further we have to show 
that E g h , s b is the steady state of the variational dynamics 
for the bath degrees of freedom, or in other words that 
the auxiliary RLM thermalizes to the groundstate with 
hybridization \fZlV k . 

This can be seen explicitly by looking at the dynamics 
of the bath degrees of freedom. To this extent we plug 
on the right hand side of equations (|30I31I) the equilib- 
rium ground-state values of W k and N kp that for a RLM 



with hybridization \fZlV k can be obtained by standard 
methods 7 . These read 



Wfl' = T^ e fa,0+ G* (iw) . V * ( 4( t ) 



as well as 



IUJ — £r 



iuj - e k 

. . „ ( 41 ) 

with V p = the impurity Green's function Q*(iuS) 

given by 



1 



(42) 



iw — Z+ A(iw) 
and the hybridization function 

V 2 

A(iw) = V fc 

k IUJ ~ £k 

If we plug these expressions into the variational equations 
of motion we can see they represent stationary solutions. 
We find indeed 



: .N kp = (£ p -e k ) TJ2 

iuj 

+T J2 e lw0+ g*(iu) 



+ V p g*{iu))V k 



(iu> - s p ) (iu> - e k ) 
' VpV k V k V p 



iuj — e k iuj — e 



= 
as well as 



(43) 



— ^ 1. 



SkVk 



iuj - e k 



+ V k (f(sk)-l) + 



+T e ^° + A*(tw)0*( 



(44) 



iuj - e k 

where in the last step we have used the fact that 

A*(iu))G*(iu>) = -l + iu)G*(iu). (45) 

This shows in particular that the groundstate value of 
the hybridization Ef lybl which immediately follows from 
W gs , is the steady state of the variational dynamics. 

1. Linearization around the steady state 

In this section we linearize the variational dynamics 
around the steady state. To this extent it is convenient to 
write the equations of motion in a different form. Indeed 
the average hybridization Eh is written as 

E hyb = V k (W k + h.c.) =2j2Vk W? (46) 
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which suggests to express the dynamics in terms of the 

, fep . A simple 



real and imaginary parts, Wj^' 1 and N] 1 ' 1 



algebra gives 



E^ 



T 

hp 



wl 



Ni p 



-- (e p - e k ) Nl p + V p Wl + % Wf 
(e p -e k ) N%-V p W?-V k W? 



(47) 
(48) 

(49) 
(50) 



These equations, together with the dynamics fo D and <p 



D = y/ZE hyb (t)si 



sirup 



p = U + E hyb {t) 



dD 



CO sip 



(51) 
(52) 



with \[Z = 4-^/ D(l/2 — D), represents the variational 

dynamics for the variables (d,4>,W^,W^,N]^,N^ 

that we want to linearize around the steady state val- 

ues, (£>*, ^, Wl, Wl,N^,NQ . 

It is a straightforward but tedious calculation to lin- 
earize this dynamics and obtain a system of first order 
linear differential equations. In principle, by computing 
the eigenmodcs of the associated matrix we could get in- 
formation of the stability of the steady state solution. 
However in order to make analytical progress it is use- 
ful to express the dynamics as a second order differential 
equation for half of the variables. This dynamics reads 

SD = a 5D + Vk SW^ (53) 



8N% = £ SI* S N? q + £ 9f bWf + lkp SD (55) 



where the coefficients a,r] k , 

> 7fcp 

read respectively 



A 



kq 



a = ^l{E* hyh f^fcos 2 ^ 



COS 

)D 

dvz: 



Vk = 2 V / Z± —qjj- cos^ip* V k E* hyh 



(56) 
(57) 



Pk = e k -Qjy cosip* V p N kp* + 



z * d -^w) E (n 2 w£ + v q v k UfflSS) 



"■'•/' = (-a ( ^ (v P w£ + v k w? 



p-k, 



z * %?r E ( w N %* + v « Vk ^*)(59) 

(60) 



Z± COStfi* V q £ q 

S kq (e 2 q + Z± cos 2 Lp± V q ) + 
- Z+cos 2 ^V q V k (61) 
n k q p q , = -{e q -e ql f 5 qp 5 q , k + (62) 

+ Z* COS 2 tp± (VpVgSqik + V k V q 5 q > p ) (63) 

= y/zlcoatpi, {S kp V p {2e q - e p ) + 6 pq V k e k )(U) 
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